Two-Dimensional Propagation of a Photoinduced Spin Wave Packet 
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We report the two-dimensional propagation of photoinduced spin wave packets in Bi-doped rare-earth iron 
garnet. Spin waves were excited nonthermally and impulsively by a circularly polarized light pulse via the 
inverse Faraday effect. Space- and time-resolved spin waves were detected with a magneto-optical pump- 
probe technique. We investigated propagation in two directions, parallel and perpendicular to the magnetic 
field. Backward volume magnetostatic waves (BVMSWs) were detected in both directions. The frequency 
of BVMSWs depends on the propagation direction. The experimental results agreed well with the dispersion 
relation of BVMSWs. 

PACS numbers: 75.50.Ee, 78.47.D-, 75.30.Ds, 78.20.Ls 



Magnetization dynamics has attracted a great deal of inter- 
est in recent years. In particular, the propagation character- 
istics of spin waves have been extensively studied because of 
their importance as the basis of magnonics, which aims to use 
packets of spin waves instead of electric currents as the infor- 
mation transmitter iulil]. Magnonics shows great promise to 
improve solid-state devices in combination with conventional 
electronics |3]. When spins are coupled mainly by the mag- 
netic dipolar interaction rather than the exchange interaction, 
the spin waves are called magnetostatic waves (MSWs) j^. 
Traditionally, MSWs have been excited by a microwave field 
emitted from a microstrip antenna. Time-resolved MSWs 
have been examined in many ways, for example, via the in- 
ductive method |5, 6], the magneto-optical Kerr effect iT Hllll . 
and Brillouin light scattering |1 2l 1 1 311 . Their nonlinear effects 
1 141] and mode quantization 1 15|] have also been investigated. 

Analysis of spin waves excited by an ultrashort light pulse 
has already been reported |'l60. In that study, the spin 
waves were excited by a heat-induced change in the sam- 
ple's anisotropy. Nonthermal excitation is thus desirable from 
the viewpoint of exciting coherent spin waves and simplify- 
ing the analysis. Recently, a novel method of magnetization 
control has been reported. The method employs the inverse 
Faraday effect (IFE) jlT'], where a circularly polarized light 
pulse nonthermally generates an effective magnetic field. IFE 
has been described as impulsive stimulated Raman scatter- 
ing BlSll . Thus, ultrafast magnetic switching is expected via 
IFE. It has been reported that IFE occurs in diverse magnetic 
materials il7l42lll . However, in these reports, discussion has 
been limited to local spin oscillations only. Non-local spin 
dynamics, i.e., propagating spin waves, excited by this im- 
pulsive excitation method has never been observed. By using 
this method, two-dimensional propagation can be observed in 
a non-contact manner, and broadband excitation can be real- 
ized, unlike resonant excitation with a microwave field. 

Yttrium iron garnet (YIG) has often been used as a sample 



in magnonics because of its intrinsically low magnetic damp- 
ing i^. In an experiment involving photoinduced spin waves, 
Bi-doped rare-earth iron garnet was found to be more suit- 
able than YIG due to its high magneto-optical susceptibility 
I22II . By doping the sample with bismuth and magnetic rare- 
earth ions, the Gilbert damping coefficient increases. Even so, 
MSWs can propagate over distances of hundreds of microme- 
ters, and they can be spatially resolved. 

Here we report the two-dimensional propagation of a spin 
wave packet excited via IFE in Bi-doped rare-earth iron gar- 
net. Propagation in two directions, parallel and perpendicular 
to the magnetic field, was investigated by an all-optical pump- 
probe experiment. 

The composition of our sample was Gd4/3Yb2/3BiFe50i2, 
and the thickness was 110 jim. The unit cell is cubic with 
three types of interstices, 16a, 24d, and 24c sites [23J . Fe^^ 
ions fill the 16a and 24d sites, and rare-earth ions fill the 24c 
sites. The magnetic moments of Fe^+ ions in different sites 
are antiferromagnetically coupled. The magnetic moments of 
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FIG. 1: (Color online) (a) The experimental configuration. The sam- 
ple surface was in the x-y plane. The magnetic field was applied 
in the x direction, (b) Time-resolved magneto-optical Faraday rota- 
tion waveforms when pump spot and probe spot overlapped spatially. 
ct(+,— ) represent pump helicities. Solid lines represent fitting re- 
sults. 
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FIG. 2: (Color online) The BVMSW band and MSSW band versus 
in-plane propagation angle 6 when Ho = 1000 Oe and //„ = 450 
Oe. The lower limit of the BVMSW band and the upper limit of 
the MSSW band moved up as Q increased. MSSW emerged when 

> Qr- 



FIG. 3: (Color online) A representative spatial map of MSWs when 
scanning parallel to the magnetic field. The delay was fixed to 2.1 ns. 
The sum of the spot diameters of the pump and probe corresponds to 
the width of the shaded region. 



rare-earth ions have much weaker antiferromagnetic coupHng 
with the net Fe^+ moment. Thus, saturation magnetization 
(\%Ms = 11 36 G) is reduced compared with the typical value 
for YIG {\nM, = 1956 G) [3]. 

Our sample had uniaxial anisotropy. The z axis was the easy 
axis, and the x-y plane was the hard plane (See Fig. 1(a)). To 
evaluate the uniaxial anisotropy field, the uniform precession 
frequency /i was determined to be 3.6 GHz under an in-plane 
external field of 1000 Oe in a magnetic resonance experiment. 
The value was compared by using 



/i = YVHo{Ho-Hu + 4nM^s), 



(1) 



where Hq is the applied magnetic field, and //„ is the uniax- 
ial anisotropy field to the z axis. Here, AnM^ff is defined as 
{N^ — Nx) 4-7tMs, Nu,\ is the demagnetizing factor along the 
x{z) axis, and 7 = 2.8 MHz/Oe is the electron gyromagnetic 
ratio. Then, //„ was determined to be 450 Oe. From the 
FWHM of the resonance peak, we also determined the Gilbert 
damping coefficient to be a = 0.02, which is about 10^ times 
higher than the typical value for YIG {^ 10^"*). This differ- 
ence was attributed to an increase in indirect coupling between 
the Fe^+ sublattice and the phonon lattice via rare-earth ions 

We employed the magneto-optical Faraday effect to mea- 
sure space- and time-resolved MSWs. The experimental con- 
figuration is shown in Fig. 1(a). Linearly polarized pulses 
from a Ti: sapphire laser with a wavelength of 792 nm and a 
pulse width of 120 fs were used as the probe. Circularly po- 
larized pulses with a wavelength of 1400 nm, generated by 
an optical parametric amplifier, were used as the pump. The 
pump beam was vertically incident on the surface of the sam- 
ple, whereas the probe beam was incident at 7°. The pump and 
probe beams were focused on the sample to spot diameters of 
about 50 and 40 /im, respectively. The spot sizes were suffi- 
ciently small to resolve MSWs spatially, as will be confirmed 
below. The photon energy of the pump pulses (0.89 eV) was 



below the lowest d-d transition in Fe^ ions (1.4 eV) II24I1 



leading to a virtual excitation. The pump fluence was 400 
mJ/cm^. In our sample, the refractive index was 2.4, and the 



absorption coefficient was 0.3 cm^' for the pump wavelength. 
Therefore, absorbed pump fluence was 1 mJ/cm'^, which en- 
abled almost uniform excitation across the entire thickness. 
An in-plane magnetic field of 1000 Oe was applied to satu- 
rate the magnetization in the x direction. The pump spot was 
scanned in the directions parallel {x axis) and perpendicular {y 
axis) to the magnetic field at intervals of 20 /im, whereas the 
position of the probe spot was fixed. Polarization rotation of 
the probe light was recorded for each pump position. When 
the pump spot was scanned in the x direction, we call it the 
parallel geometry. When the pump spot was scanned in the y 
direction, we call it the perpendicular geometry. All measure- 
ments were performed at room temperature, which is below 
the Curie temperature of our sample (7^ = 573 K). 

Figure 1(b) shows polarization rotation versus the time de- 
lay between the probe pulse and the pump pulse when the 
two spots overlapped spatially. The phase of spin oscillation 
shifted by n rad when the helicity of the pump light reversed. 
This is clear evidence of IFF, where the direction of the mag- 
netic field is determined by the pump helicity. The damped 
oscillations were fitted by Aexp(— f/T)sin(27i:/f) B. The 
fitting results yielded A == 5.3 mrad, B = 1.8 mrad, T = 0.52 
ns, and f — 2.9 GHz. Because the direction of the effec- 
tive magnetic field pulse generated via IFF was perpendicu- 
lar to the sample surface, the oscillations showed sinusoidal- 
dependence. In a sample having a Gilbert damping coeffi- 
cient a, the spin precession is damped in the direction of 
the magnetic field on a time scale of l/aco [4J- In our case, 
l/aco = 2.7 ns is considerably longer than the value of T de- 
termined by the fitting. This suggests that the energy and an- 
gular momentum propagated away from the pump spot in the 
form of spin waves. The relatively long-tailed (over 3 ns) de- 
crease of polarization rotation, represented by B, was indepen- 
dent of the pump helicity and was attributed to a heat-induced 
change in anisotropy. Nevertheless, IFF clearly dominated the 
excitation of spin oscillations within 3 ns. This thermal effect 
is subtracted in later analyses. Hereafter, the pump helicity 
was fixed to C7(+). 

In an in-plane magnetized plate, two types of MSWs ex- 
ist: backward volume magnetostatic waves (BVMSWs) and 
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FIG. 4: (Color online) Spin waves when scanning parallel to the mag- 
netic field, (a) A spatiotemporal map of spin waves. The color bar 
represents polarization rotation of the transmitted probe light in units 
of mrad. The dotted line represents calculation results of group ve- 
locity based on Eq. (4). (b) Dependence of Fourier spectrum of 
MSWs on scan length. The dotted lines respectively represent calcu- 
lated upper and lower limits of the B VMSW band. 



magnetostatic surface waves (MSSWs) 112611 . Hereafter, we 
discuss MSWs in a uniaxial anisotropic plate with perpendic- 
ular easy axis. In-plane propagation angle 9 is defined as the 
angle between the direction of the magnetic field (x axis) and 
the in-plane spin-wavevector BVMSWs propagate at any an- 
gle to the magnetic field, whereas MSSWs propagate at an- 
gles greater than a critical angle 9c = sin" ' ( ^/tHq/Ji 1 to the 



magnetic field |27]. These MSWs exist in a finite range of fre- 
quencies, which are called the BVMSW band and the MSSW 
band, respectively. Figure 2 shows the frequency range of the 
bands versus in-plane propagation angle 9 for Ho — 1000 Oe 
and //„ = 450 Oe. When 9 = 0°, the BVMSW band exists in 
a range /o < / < /i ■ Here, /o is defined as 



/o = rVHo{Ho-H,). 



(2) 



The lower limit of the BVMSW band moves up from /o to 
Yy/{Ho — Hu) (//() + 4;rMeff) as the propagation angle 9 in- 
creases from 0° to 90°. MSSWs emerge when > 0,. = 51°, 
and the upper limit of the MSSW band moves up from /i to 
/2 as the propagation angle 9 increases from 9c to 90°. Here, 
/2 is defined as 



fi 



/ y(47rMeff-//„) 




(3) 



It should be noted that spin waves within a finite range of 
propagation angles 9 were detected due to the finite sizes of 
the pump and probe spots in our experiment. Therefore, the 
detected range of 9 depended on the scan length, which is de- 
fined as the distance between the pump and probe spots. This 
effect caused the spectral shape of MSWs to be dependent on 
the scan length, as will be mentioned below. 

Figure 3 is a representative spatial map of spin waves in the 
parallel geometry. The delay time was fixed to 2.1 ns. In this 
geometry, propagation of BVMSWs is expected ||4l|26|]. In the 
shaded region, whose width corresponds to the sum of the spot 
diameters of the pump and probe, the thermal effect cannot 



FIG. 5: (Color online) Spin waves when scanning perpendicular to 
the magnetic field, (a) A spatiotemporal map of spin waves. The 
color bar represents polarization rotation of the transmitted probe 
light in units of mrad. (b) Dependence of Fourier spectrum of MSWs 
on scan length. The dotted lines respectively represent calculated up- 
per and lower limits of the BVMSW band and the MSSW band. 



be completely subtracted. Outside this region, propagation of 
spin waves with a wavelength A = 300 /im was observed. The 
spot sizes were confirmed to be small enough to resolve the 
waveform. The wavenumber (k = 2n/A,) also satisfies ko ^ 
k <^ l//o. Here, ko is the wavenumber of the electromagnetic 
waves that possess the same frequency as the spin waves, and 
lo — ^Dj'XnMs is an exchange length. Using a value of the 
exchange stiffness from the Uterature, D = 9.4 x 10"^ Oe cm^ 
ll28ll . /q is 29 nm. The exchange interaction is negligible when 
the above condition is satisfied. Therefore, these waves are 
regarded as being MSWs. The propagation in both directions 
was symmetric in amplitude and phase, unlike the excitation 
by a microwave field emitted from a microstrip antenna 12911 . 
This is because the direction of the effective magnetic field by 
IFF is perpendicular to the sample surface, which affects the 
magnetization symmetrically about the pump spot. 

Figure 4(a) shows a spatiotemporal map of the spin waves 
in the parallel geometry. It is clearly seen that a spin wave 
packet propagated from the pump spot, and the phase of the 
wave packet moved toward the pump position from the out- 
side, which is a characteristic of BVMSWs. The dispersion of 
the lowest-order mode of the BVMSWs is approximated by 



a 



/w 



Ho{Ho-H^,+4-7zMefi 



1 — exp(— fed) 
kd 



(4) 



Here, li = 110 fim is the thickness of the sample. When we use 
the experimentally determined wavenumber (2;r/300 /xm"') 
of the MSWs, / is calculated to be 2.8 GHz from Eq. (4), 
which is consistent with the experimentally determined fre- 
quency (2.9 GHz). The group velocity of the BVMSWs cal- 
culated from Eq. (4) is 135 km/s, which is shown as the dotted 
line in Fig. 4(a). This result is in good agreement with the ex- 
perimental results. Figure 4(b) shows the dependence of the 
Fourier spectrum of the MSWs on the scan length in the par- 
allel geometry. The upper limit /i and the lower limit /o of 
the BVMSW band, which are shown in Fig. 4(b), were cal- 
culated from Eqs. (1) and (2). The frequency range of the 



propagating MSWs, centered on 2.9 GHz, agreed well with 
the calculated BVMSW band (see Fig. 2). As the pump spot 
went away from the probe spot, only the propagation angle 9 
close to 0° was detected. The BVMSW band was the broadest 
for = 0°, in agreement with Fig. 2. 

Figure 5(a) shows a spatiotemporal map of the MSWs in 
the perpendicular geometry, where propagation of MSSWs 
was also expected in addition to BVMSWs. However, the 
phase of the wave packet went toward the pump position from 
the outside. This suggests that the BVMSWs were the dom- 
inant MSWs propagating in this direction. As the pump spot 
went away from the probe spot, only the propagation an- 
gle close to 90° was detected in this geometry. Because 
the BVMSW band width was reduced as the scan length in- 
creased, the group velocity of the BVMSWs became small in 
this direction. Unlike the parallel geometry, the wave packet 
stayed near the pump spot due to this effect. Phase discon- 
tinuity around the scan length of 40 fim is an artifact which 
is attributed to partial spatial overlap between the pump and 
probe spots. Figure 5(b) shows the dependence of the Fourier 
spectrum of the MSWs on the scan length in the perpendic- 
ular geometry. The spectral intensity at the lower frequency 
in the BVMSW band was suppressed more than that at the 
higher frequency as the scan length increased, because the 
lower limit of the BVMSW band moved up as the scan length 
increased. The upper limit of the MSS W band, /2, in Fig. 5(b) 
was calculated from Eq. (3). The amplitude of the MSSWs 
was very small compared with that of the BVMSWs, which 
suggests that the MSSWs were not effectively excited by the 
pump pulses. Lenk et al. |16|] suggested that asymmetric ex- 
citation across the thickness of the sample is needed to excite 
MSSWs. Because of the neghgible absorbance of the pump 
pulse, asymmetric excitation was not realized in our sample. 
Thus, the excitation efficiency of MSSWs was lower than that 
of BVMSWs, which were distributed almost uniformly across 
the thickness of the sample 02711 . 

In summary, two-dimensional propagation of photo- 
induced spin waves in Bi-doped rare-earth kon garnet was 
investigated by an all-optical pump-probe experiment. Prop- 
agating spin waves were impulsively excited by a circularly 
polarized light pulse via IFF. In the parallel geometry, BVM- 
SWs were detected, and their frequency and wavelength were 
in good agreement with the theoretical dispersion relation. 
BVMSWs were also detected in the perpendicular geome- 
try, with a higher frequency than that in the parallel geome- 
try. This frequency dependence on the propagating direction 
is consistent with the theoretical prediction that the lower limit 
of the BVMSW band moves up as the propagation angle in- 
creases. MSSWs were not effectively excited, which can be 
attributed to the negligible absorbance of the pump pulse in 



our sample. Our findings open up the possibility of using IFF 
as an impulsive spin wave generator. This nonthermal, coher- 
ent excitation method will be an important tool in magnonics. 
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